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We present a core-collapse supernova model for the extremely luminous Type Ic supernova 2007bi. 
By performing numerical calculations of hydrodynamics, nucleosynthesis, and radiation transport, we 
find that SN 2007bi is consistent with the core-collapse supernova explosion of a 43 Mq carbon and 
oxygen core obtained from the evolution of a progenitor star with a main sequence mass of 100 Mq and 
metallicity of ^ = Zq/200, from which its hydrogen and helium envelopes are artificially stripped. 
The ejecta mass and the ejecta kinetic energy of the models are 40 Mq and 3.6 x 10^^ erg. The ejected 
^^Ni mass is as large as 6.1 Mq, which results from the explosive nucleosynthesis with large explosion 
energy. We also confirm that SN 2007bi is consistent with a pair-instability supernova model as has 
recently been claimed. We show that the earlier light curve data can discriminate between the models 
for such luminous supernovae. 

Subject headings: supernovae: general, supernovae: individual (SN 2007bi, SN 2006gy), gamma-ray 
burst: general 



1. INTRODUCTION 

A massive star with the main sequence mass (Mms) 
in the range of 10 — 140 Mq forms an Fe core in its 
center and eventually collapses. This collapse is thought 
to end up with the core-collapse supernova (SN) of Type 
II, lb, or Ic (Filippenko 1997 for a review). If a star 
is as massive as Mms= 140 — 300 Mq, the oxygen-rich 
core becomes dynamically unstable owing to the electron- 
positron pair creations (Rakavy & Shaviv 1967; Barkat, 
Rakavy, & Sack 1967). As the internal energy is spent 
by the pair creations, the core loses the stability and 
starts to collapse. When the central temperature exceeds 
^ 5 X 10^ K, the core becomes stable but the temper- 
ature is so high that oxygen burning becomes explosive 
to produce enough energy to unbind the star. A large 
amount of ^^Ni is synthesized by the explosive burning 
{e.g., Umeda & Nomoto 2002, hereafter UN02; Heger 
& Woosley 2002) and the subsequent radioactive decays 
power the light curve. Thus this event is theoretically 
predicted to be observed as a pair-instability supernova 
(PISN). Some luminous SNe like SN 2006gy (Ofek et 
al. 2007; Smith et al. 2007) have been suggested to be 
PISNe (see Section 4.2), but no clear consensus has been 
reached {e.g., Kawabata et al. 2009). 

Recently, Gal- Yam et al. (2009) (G09 hereafter) sug- 
gested that the extremely luminous Type Ic SN 2007bi 
is the first observed example of the PISN. They showed 
that the PISN model is consistent with the observed 
light curve (LC) and the nebular spectra of SN 2007bi. 
They estimated the masses of C, O, Na, Mg, Ca, and 
^^Ni from the observed optical spectra. Other elements 
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with no strong emission lines in the optical range. Si 
and S, are assumed to be the same as the PISN model of 
Heger & Woosley (2002). Young et al. (2010) (YIO here- 
after) showed multi-color observations of SN 2007bi and 
the metallicity of the host galaxy. However, the above 
observations of SN 2007bi have not quantitatively been 
compared with the core-collapse SN models. In view of 
the importance of clarifying the final fates of very mas- 
sive stars, we examine how strong the observational con- 
straint on the theoretical models are. 

The aim of this Letter is to show that a core-collapse 
SN model is indeed consistent with the observed prop- 
erties of SN 2007bi, if the progenitor is as massive as 
Mms~ 100 Mq and the explosion energy is large. This 
would imply that SN 2007bi might not necessarily be a 
PISN. In Section 2, we summarize the progenitor models 
for SN 2007bi and numerical methods used for our calcu- 
lations of hydrodynamics, nucleosynthesis, and the LC. 
The core-collapse SN models of SN 2007bi are presented 
in Section 3 and the results are discussed in Section 4. 

2. PROGENITOR AND EXPLOSION MODELING 
2.1. Progenitor 

The high peak luminosity and the long rise time of the 
LC of SN 2007bi (G09, YIO) require a large amount of 
^®Ni (> 3 Mq, G09) and a large ejecta mass. These ob- 
servations imply that the progenitor of SN 2007bi is mas- 
sive. We apply a pre-SN model with Mms = 100 Mq cal- 
culated by Umeda & Nomoto (2008, UN08 hereafter). 
UN08 assumed the metallicity of the progenitor models 
to be Z = Zo/200, which is small enough to avoid a 
large amount of wind mass loss. Then the pre-SN model 
remains as massive as M = 83 Mq, whose carbon -I- oxy- 
gen (C-l-0) core is massive enough (43 Mq) to produce 
a large amount of ^^Ni. 

However, the pre-SN model has a massive H-rich enve- 
lope, while SN 2007bi does not show the lines of either 
H or He. Therefore the progenitor must have lost its H- 
rich envelope (36 Mq) and He layer (4 Mq) during the 
pre-SN evolution, thus having only the bare C-l-0 core at 
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Fig. 1.— Bolometric LCs of the C+O star SN models CClOO 
(Mej = 40 Mq, Ekin = 3.6 X 10^2 gj,g^ and MseNi = 6-1 Mq). 
The observed bolometric LC (open circles) is taken from YIO. The 
bolometric magnitude of the rise part of SN 2007bi (open square) is 
estimated from the R band magnitude. All the calculated LC has 
the same physical structure but the degrees of mixing are different. 
The horizontal axis shows the days in the rest frame. 



the explosion. We construct the pre-SN C+O star model 
of 43 Mq, by removing the H-rich envelope and He layer 
from the 83 Mq star. Note that the metallicity of the 
host galaxy of SN 2007bi {Z - Zq/3, YIO) is higher 
than that of our adopted progenitor [Z = Zq/200). The 
wind mass loss is expected to work more efficiently and 
the main sequence mass of the progenitor which has the 
C+O core mass of 43 Mq might be more massive. The 
rotation of stars can also play a role in the mass loss 
[e.g., Meynet et al. 2003; Hirschi et al. 2004; Maeder et 
al. 2005; Georgy et al. 2009). Another possible cause 
of such envelope stripping is the formation of a com- 
mon envelope during a close binary system, where the 
smaller mass companion star spirals into the envelope of 
the more massive star. The outcome depends on whether 
the energy available from the spiral-in exceeds the bind- 
ing energy of the common envelope, thus being either a 
merging of the two stars or the formation of two compact 
stars, e.g., a C+O star and a He star. 

2.2. Explosion Modeling 

We calculate the explosion of the pre-SN C+O star 
(43 M0) as described above. Explosions are induced 
by a thermal bomb and followed by a one-dimensional 
Lagrangian code with the piecewise parabolic method 
(Colella & Woodward 1984). Note that the explosion 
energy is a free parameter in core-collapse SN explosion 
models while it is not in PISN explosion models. Explo- 
sive nucleosynthesis is calculated as post-processing for 
the thermodynamical history obtained by the hydrody- 
namical calculations. The resultant abundance distribu- 
tion is basically very similar to those calculated by UN08 
(see Figures 5 and 6 in UN08). 

The dynamics of the ejecta is followed until 1 day af- 
ter the explosion, when the expansion already becomes 
homologous (r (X w). The bolometric LCs are calculated 
for the homologous ejecta by using the LTE radiation 
transfer code (Iwamato et al. 2000) that includes the ra- 
dioactive decays of ^®Ni and ^^Co as energy sources. This 
code calculates the 7-ray transport for a constant 7-ray 
opacity (0.027 cm^ g~i, Axelrod 1980) and assumes all 
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Fig. 2. — Photospheric velocities of the core-collapse SN and 
PISN models. The horizontal axis shows the days in the rest frame. 
The line velocities of SN 2007bi observed by YIO are also shown for 
comparison. The line velocities of O I A7774 shows the lower limit. 
All the photospheric velocities are low enough to be consistent with 
the observed line velocities. 

the emitted positrons are absorbed in situ^ . For the op- 
tical radiation transport, the Thomson scattering opac- 
ity is obtained by calculating the electron density from 
the Saha equation, and the Rosseland mean opacity is 
estimated from the empirical relation to the Thomson 
scattering opacity (Deng et al. 2005). 

3. CORE-COLLAPSE SUPERNOVA MODELS FOR SN 2007BI 

We construct several core-collapse SN models and com- 
pare with the observations of the bolometric LC and the 
line velocities of SN 2007bi shown in YIO. Since the LC 
of YIO does not cover the rising part of the LC, we es- 
timate the bolometric magnitude of the rising part from 
the R band observations (G09) assuming the same bolo- 
metric correction (0.45 mag) as in the R band maximum. 
We also take into account mixing since it is possible that 
a jet emerges from the central remnant and causes the 
mixing of the ejecta [e.g., Maeda & Nomoto 2003, Tom- 
inaga 2009). 

The LCs of successful models (CClOO) are shown in 
Figured] The kinetic energy (Skin), ejecta mass (Afoj), 
and ^^Ni mass [M&e^-^) in the ejecta are -Ekin = 3.6 x 10^^ 
erg, Mej = 40 Mq, and Msgn; = 6.1 Mq, which are the 
same in all the models. The mass cut between the ejecta 
and the compact remnant is set at Mr — 3 -Mq, where 
Mr is the mass coordinate, so that the ejecta contains 
6.1 Mq of ^^Ni, which turns out to be consistent with 
the bolometric LC of SN 2007bi. The kinetic energy 
-Ekin needs to be large to produce MseNj = 6.1 Mq. The 
mass of some elements in the ejecta are summarized in 
Table [T] The kinetic energy is as large as those of previ- 
ously observed SNe which were associated with a gamma- 
ray iDurst (SNe 1998bw, 2003dh, 20031w; e.g., Nomoto et 
al. 2006). 

We adopt two different degrees of mixing to see its 
effects on the LC. The full- mixing model assumes that 
the whole ejecta are uniformly mixed. The half-mixing 
model assumes that the inner half of the ejecta (in the 
mass coordinate) is uniformly mixed. One of the effects 
of the mixing is seen in the rise time of the LC. With 

^ This assumption of the positron absorption does not have much 
effect on the LCs we show in this paper, because the contribution 
from the gamma-rays is still a dominant energy source of them. 
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TABLE 1 

The Amount of Elements Contained in the Ejecta 



12c 20Ne 24]v[g 


28 Si 


32g 


36Ar 40ca 56Ni 


1.4 18.7 1.4 1.5 


5.1 


2.7 


0.5 0.4 6.1 



Note. — Units: Mq 



mixing, ^^Ni is distributed closer to the surface of the 
ejecta, so that the diffusion time is shorter and the rise 
time becomes shorter. The rise time of the model with- 
out mixing is 85 days while the rise times of the half- 
mixing model and the full-mixing model are 67 days and 
52 days, respectively. As the rise time of SN 2007bi is 
not observationally well-determined, all the models are 
consistent with the bolometric LC of SN 2007bi. The ini- 
tial decline part of the calculated LCs before maximum 
is formed by the shock heating of the envelope and its 
subsequent cooling due to rapid expansion. Radiation 
hydrodynamical calculations are required to obtain the 
realistic LC at this epoch. 

In Figure [21 we show the photospheric velocities ob- 
tained by the LC calculations. With the photospheric 
velocities, we also show the observed line velocities of SN 
2007bi taken from Figure 17 of YIO. The photospheric 
velocities of all the models are consistent with the ob- 
served lowest line velocities, which are thought to trace 
the photospheric velocity. 

One of the big difference between the core-collapse SN 
models and the PISN models is the abundance of the el- 
ements like Si and S. The abundance of our core-collapse 
SN model is consistent with the directly estimated abun- 
dances from the one-zone model of G09 adopted to the 
observed emission lines, i.e., C, O, Na, Mg, Ca, and 
^^Ni. However, it differs significantly in elements whose 
abundances are only indirectly constrained by the spec- 
tra (such as Si and S) as these may have played a role in 
line cooling processes. In order to confirm that the abun- 
dance of the core-collapse SN model are consistent with 
the nebular spectra, we have to perform spectral synthe- 
sis calculations for the realistic hydrodynamical model of 
ejecta rather than the single-zone adopted by G09. As Si 
and S have many emission lines in the infrared range, in- 
frared spectra are also helpful to distinguish PISNe from 
core-collapse SNe. We also point out that, if SN 2007bi 
is confirmed to be a PISN, we could expect that PISNe 
played a role in the chemical enrichment in the early Uni- 
verse and there should be some old stars with chemical 
compositions expected from PISNe, although they are 
still not discovered {e.g., Cayrel et al. 2004). 

4. CONCLUSIONS AND DISCUSSION 

In this Letter, we have shown that the LC and pho- 
tospheric velocity of SN 2007bi are well-reproduced by 
the core-collapse SN model CCIOO. As some gamma-ray 
bursts are connected to such high energy Type Ic SNe, 
the extremely luminous SNe like SN 2007bi could also be 
connected to gamma-ray bursts which result from very 
massive stars. If this is the case, extremely luminous SNe 
like SN 2007bi could be connected to gamma-ray bursts 
of much more massive star origin than known SNe asso- 
ciated with a gamma-ray burst. Even stars more massive 
than 300 Af© could be the origin of luminous SNe {e.g., 
Ohkubo et al. 2006, 2009). 
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Fig. 3.— A PISN model for SN 2006gy {PISN270). The bolo- 
metric LC of SN 20D6gy is taken from YIO and Agnoletto et al. 
(2009). The first point of SN 2006gy is the R band upper limit 
from Smith et al. (2007). The CCIOO no-mixing model and ob- 
served bolometric LC of SN 2007bi in Figure[l]are also shown. The 
horizontal axis shows the days in the rest frame. 

We note, however, that, although SN 2007bi may not 
necessarily be a PISN, the observational data available 
for SN 2007bi is not sufficient to single out the ex- 
plosion mechanism. In fact, Kasen & Bildsten (2009) 
suggested that the magnetar-powered light curve model 
(also, Maeda et al. 2007; Woosley 2009) might explain 
the LC of SN 2007bi. 

Here we show the comparison between our PISN model 
and SN 2007bi, and discuss how to distinguish the models 
for luminous SNe. We also apply such LC comparison to 
SN 2006gy. 

4.1. Pair-Instability Supernova Models for SN 2007bi 

In Section 3 (Figure [Ij, we have shown that obser- 
vations of SN 2007bi are well reproduced by the core- 
collapse SN model (CCIOO). Here we confirm the claim 
made by G09 that a PISN model can also be consistent 
with the bolometric LC of SN 2007bi by using the ap- 
proximate PISN model PISN270. 

The PISN270 model is constructed by scaling the 
physical structure of the homologously expanding model 
CCIOO to the ejecta model with M^j = 121 Mq and 
Ekin — 7 X 10^^ erg. The ejecta mass Moj is the 
same as the C-l-0 core mass of the PISN model with 
Mms = 270 Mq (UN02), and E-i.in is obtained from 
the nuclear energy released by explosive nuclear burn- 
ing of the C-l-0 core (UN02). Here the same amount 
of ^^Ni (M56Ni = 9.8 Mq) as in the 270 Mq model 
(UN02) is assumed to be synthesized in the inner layers. 
Note that the 270 Mq model of UN02 still has the H- 
rich and He envelopes at the time of explosion and, here, 
we assume that the envelopes were stripped off by some 
mechanism. 

Figure [3] shows that the bolometric LC of PISN270 
(the red line) is consistent with the bolometric LC of 
SN 2007bi (red open circles). The rise time to the LC 
peak for PISN270 is ^ 150 days, being consistent with 
the PISN model in G09. This rise time is longer than 
the core-collapse SN model CCIOO (Figure [3]), because 
the photon diffusion takes more time in more massive 
PISN270. Ahhough MseNi of PISN270 is - 1.6 times 
larger than that of the core-collapse SN CCIOO model, 
the longer rise time lowers the peak brightness powered 
by the radioactive decay. These two effects make the 



4 



Moriya et al. 



peak magnitude of PISN270 similar to CCIOO. 

This difference in the rising part of the LC is impor- 
tant to discriminate between the core-collapse SN and 
the PISN models. Although SN 2007bi was not observed 
early enough, much earlier observations before the peak 
could constrain the SN type from the LC. In addition, 
as already mentioned in Section 3, the abundance of Si 
and S would also be a key to distinguish between the two 
models. 

4.2. Models for SN 2006gy 

As mentioned in Section 1, there has been some sug- 
gestions that the luminous Type Iln SN 2006gy is a PISN 
[e.g., Smith et al. 2007). We thus apply our LC mod- 
els for comparison with SN 2006gy. Figure [3] shows the 
bolometric LC of SN 2006gy (filled circles and triangles). 
As only the R band magnitude was observed in the early 
epochs of SN 2006gy (the filled square), we cannot con- 
struct the bolometric LC at the early epochs but we can 
constrain the rise time of the LC. Our calculations show 
that the rise time of the PISN model is too slow to be 



consistent with SN 2006gy. Although our PISN model 
does not have a H-rich envelope, the presence of the H- 
rich envelope could even slow down the brightening {e.g., 
Kawabata et al. 2009). Woosley et al. (2008) showed 
that the interaction between the pulsating core and the 
envelope can power the LC of SN 2006gy. As a similar 
mechanism, the interaction of a SN ejecta with its very 
dense circumstellar matter could convert the kinetic en- 
ergy of ejecta directly to radiation energy and could also 
be the origin a luminous SN like SN 2006gy. 
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